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1. Introduction 



Using lattice QCD to compute the nonperturbative contributions to 5-physics quantities allows 
us to strengthen constraints on the CKM unitarity triangle and determine CKM matrix elements. 
In the game of over-constraining the apex of the CKM unitarity triangle, 5-physics may also reveal 
signs of physics beyond the Standard Model. The standard global fit of the CKM unitarity triangle 
uses 5-meson decay constants and mixing as an input [1-3].^ Experiments such as BaBar, Belle, 
and LHCb measure Bq — Bq mixing in terms of mass differences (oscillation frequencies) Mdq to 
subpercent accuracy [4]. In the Standard Model, these are parameterized by [5] 

AM, = 9^y^^SoMBjlBB^K^V4\ (1.1) 

where the index q denotes a J- or a ^-quark and Mb^ the mass of the 5^-meson. Perturbatively 
computed inputs are the QCD coefficient T]^ [5] and the Inami-Lim function 5o [6]; contributions to 
be computed nonperturbatively are the leptonic 5-meson decay constant fs^ and the bag parameter 
Bsq. Precise theoretical determinations of /b^Bb^ are needed for precise extraction of the product 
of CKM matrix elements VfqVtb- It is in particular advantageous to compute the ratio [7] 

because statistical and systematic uncertainties largely cancel and the nonperturbative contribution 
is contained in the SU (3) breaking ratio 

^^ Ib^^Bb^ (1.3) 

fBd\/BBd 

Further, constraints of the CKM triangle from BR{B tv) require precise knowledge of the decay 
constant [8-10]. Currently, only HPQCD and Fermilab/MILC have published results for the 
decay constants and /^^ as well as t, using 2+1 dynamical flavor gauge field configurations and 
ensembles with different lattice spacing. However, both groups use MILC's Asqtad lattices [1 1-14] 
and hence an independent cross-check is desirable. 

Simulating 5-physics on the lattice provides a special challenge because an additional scale 
given by the large Z^-quark mass needs to be accommodated. Here we report on our project to 
compute 5-physics quantities using the nonperturbatively tuned relativistic heavy quark action for 
Z7-quarks and domain-wall fermions for the light w, J, ^-quarks. 



2. Computational setup 

We base our 5-physics project on gauge-field configurations generated by the RBC/UKQCD 
collaborations with 2+1 flavors of domain- wall fermions and the Iwasaki gluon action [15, 16]. The 
ensembles and their parameters are listed in Tab. 1. Using the domain-wall action [17, 18] we gen- 
erate on each configuration six light valence-quark propagators with quark masses ^^m^^j = 0.005, 

^ For the latest results see: http://ckmfitter.in2p3.fr/, http://utfit.romal.infn.it/ and 
www . lattice saver ages . org. 
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Table 1: Lattice simulation parameters used in our 5-physics program. The columns list the lattice volume, 
approximate lattice spacing, light (m/) and strange (m/^) sea-quark masses, unitary pion mass, and number of 
configurations and time sources analyzed. 



# time 



{L/af > 


; (T/a) 


^ a(fm) 


ami 


anth 


M;,(MeV) 


# configs. 


sources 


) 


<64 


0.11 


0.005 


0.040 


329 


1636 


1 


243 ) 


<64 


0.11 


0.010 


0.040 


422 


1419 


1 


32^ ) 


<64 


0.086 


0.004 


0.030 


289 


628 


2 


32^ ) 


<64 


0.086 


0.006 


0.030 


345 


889 


2 


32^ ) 


<64 


0.086 


0.008 


0.030 


394 


544 


2 



0.010, 0.020, 0.030, 0.0343 and 0.040 on the coarser 24^ ensembles and amf^^ = 0.004, 0.006, 
0.008, 0.025, 0.0272 and 0.030 on the finer 32^ ensembles. On the two "coarse" 24^ ensembles 
(a^O.l Ifm; = 1 .729 GeV) with light sea-quark mass ami = 0.005 and 0.01 we place one time 
source per configuration, while on the "fine" 32^ ensembles (a ^ 0.086 fm; = 2.281 GeV) with 
light sea-quark masses ami = 0.004, 0.006 and 0.008 we place two time sources per configuration 
separated by half the temporal extent of the lattice. The mass of the light quark propagators with 
^^vai ^ 0.0343 and am^^ = 0.0272 is extremely close to the value of the physical strange quark 
mass on the corresponding ensembles: am^^ = 0.0348(11) and am^^ = 0.0273(11). 

For the b quarks we use the anisotropic Sheikholeslami-Wohlert (clover) action with the rela- 
tivistic heavy-quark (RHQ) interpretation [19, 20]. We require two experimental inputs for tuning 
the three parameters, moa, cp, ^, nonperturbatively. In our notation the RHQ action is defined by 



mo + 70^0 ^ + C 



a{D^^ 



y-D- 



ICp 



I n 



(2.1) 



We tune the parameters by probing seven points of the parameter space (see Fig. 1) measuring each 
time the spin averaged 5^-meson mass M, the hyperfine splitting Am and the ratio of the rest mass 
over the kinetic mass M1/M2. Assuming linear dependence of these measured quantities on the 
parameters for the seven points used, we obtain the tuned parameters by matching M and Am to 
experimental results and demanding that M1/M2 equals one. Typically this requires a few iterations 
to ensure that the tuned parameters are given by an interpolation in a region where the measured 
parameter dependence is sufficiently linear. Further details are presented in [21] where we also 
predict bottomonium masses and mass splittings as test of our tuned parameters. 

The computation of the 5-meson decay constants presented here is performed similarly to the 
tuning, i.e. we compute decay constants at the same seven points in parameter space (see Fig. 1) 
and then interpolate these measurements to the value of the tuned RHQ parameters given in Tab. 2. 
We measure the decay constant /b by computing on the lattice the decay amplitude which is 
proportional to the vacuum-to-meson matrix element of the heavy-light axial vector current ^ = 
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Figure 1: The seven RHQ parameter sets are indicated by the six endpoints and the central point of the 
molecule-like structure in the m^a, cp and ^ parameter space. 



Table 2: Tuned parameters m^a, cp, ^ of the RHQ action for simulating /?-quarks on the ensembles listed in 
Table 1. For details see [21]. 



f» a{fm) 


'"'sea 


niQa 


Cp 


c 


0.11 
0.11 


0.005 
0.010 


8.43(7) 
8.47(9) 


5.7(2) 
5.8(2) 


3.11(9) 
3.1(2) 


0.086 
0.086 
0.086 


0.004 
0.006 
0.008 


4.07(6) 
3.97(5) 
3.95(6) 


3.7(1) 
3.5(1) 
3.6(1) 


1.86(8) 
1.94(6) 
1.99(8) 



{Q\£^^\B{p))/^B = ip^^f/MB (2.2) 

with the four-momentum of the 5-meson. To reduce discretization errors we use an 0{a)- 
improved axial vector current on the lattice 

where the coefficient ci is computed at 1-loop with mean-field improved lattice perturbation theory 
[22]. IMultiplying the decay amplitude by the renormalization factor Z^, the lattice spacing a and 
the mass of the 5^-meson Mb^ we obtain the decay constant 

fB,^Z^^^'a-'l^/^^. (2.4) 

In our final setup we intend to compute the renormalization factor following the mostly non- 
perturbative method outlined in [23], i.e. is split into a nonperturbative part containing the 
flavor-conserving factors Zy and Z^ and a perturbatively computed factor p^/ which is expected to 
be close to one and to have a more convergent series expansion in a^: 

Z^^Pu^ZfzJl (2.5) 
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Figure 2: Renormalized and -improved decay amplitude computed on the full data set and shown 
in (232 lattice units. For better visibility data points for m^^^ = 0-01 and 0.008 are plotted with a small, 
positive, horizontal offset; data points for ml^^ = 0.004 have a small, negative, horizontal offset. The vertical 
gray lines with error-band indicate the physical values of the d- and ^-quark masses [15, 16, 25]. 



Here we use the nonperturbatively obtained Zy [16], while we compute Zy and p^i at 1-loop in 
mean-field improved lattice perturbation theory [22]. Work is in progress to obtain Zy^ nonpertur- 
batively on all ensembles [24]. 



3. Preliminary results 

In Fig. 2 we show our preliminary results for the renormalized, 0((2) -improved decay ampli- 
tude = Z(j)^^^. While a chiral extrapolation to the physical J-quark mass is needed to obtain 
<I>5^, we may take advantage of the simulation data close to the physical strange quark mass to read 
off a rough estimate of ^b^- This indicates our result is in agreement with the existing literature and 
we expect statistical uncertainties of a few MeV. Figure 3 shows our results for the SU(3) breaking 
ratio of ^Bs/^Bg- The data exhibit a weak dependence on the light sea-quark mass as well as a 
mild dependence on the lattice spacing. 



4. Outlook 

Currently work is in progress on the nonperturbative computation of Zy' as well as checks 
on the perturbative computation of the correction factor p^/ and the 0(a) -improvement coefficient 
ci. We are also woridng on combined, correlated fits using expressions based on SU(3) and SU(2) 
heavy-light meson chiral perturbation theory [26-29]. With the individual parts being in good 
shape we hope to publish results soon. We will move on to the computation of Bq — Bq mixing 
matrix elements and their ratio t, . 
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Figure 3: Ratio of decay amplitudes ^Bs/^Bq computed on the full data set and shown in a^2 lattice units. 
For better visibility data points for ml^^ = 0.01 and 0.008 are plotted with a small, positive, horizontal offset; 
data points for m^^^ = 0.004 have a small, negative, horizontal offset. The vertical gray lines with error-band 
indicate the physical values of the d- and ^--quark masses [15, 16, 25]. The open black triangle (square) 
marks the value of the valence strange-quark mass in the simulation at which the plotted ratio is one by 
construction. 
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